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Evolution of a neutrino-cooled disc in Gamma- Ray Bursts 
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ABSTRACT 

Rapid, hyper-Eddington accretion is likely to power the central engines of gamma-ray 
bursts (GRBs). In the extreme conditions of densities and temperatures the accreting 
torus is cooled by neutrino emission rather than by radiation. Another important 
cooling mechanism is the advection of energy into the central black hole. We compute 
the time evolution of a neutrino-dominated disc that proceeds during the burst and 
investigate the changes in its density and temperature. The discrimination between 
short and long bursts is made on the basis of the different rates of material inflow 
to the outer parts of the disc, thus favoring the binary merger scenario for the short 
GRBs and the coUapsar scenario for the long ones. Within the context of the coUapsar 
model, we also study the evolution of the photon luminosity of the remnant disc up to 
times ~ day, and discuss its implications for the production of emission lines in GRB 
spectra. 
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1 INTRODUCTION 

In the current view, GRBs are thought to be produced in rel- 
ativistic ejecta that dissipate energy by internal shocks (Rees 
& Meszaros 1994). However, other gamma-ray production 
mechanisms have also been proposed; these include magnetic 
field reconnections (Drenkhahn & Spruit 2002; Sikora et al. 
2003) in Poynting-fiux dominated explosions (Lyutikov & 
Blandford 2003), thermal photospheric emission (Paczyriski 
1986; Thompson 1994), Compton drag (Lazzati et al. 2000) 
or photon-pair cascades and multiple Compton scatterings 
(Stern 2003). The proposed origin of the relativistic ejecta 
invokes the birth of a black hole in a catastrophic event. 
This is thought to take place after the merger of two neu- 
tron stars or a neutron star and a black hole (Eichler et al. 
1989; Paczyriski 1991; Narayan, Paczyriski & Piran 1992) or 
in a "coUapsar" (Woosley 1993; Paczyriski 1998). Both these 
models involve a stage in which the black hole is surrounded 
by a debris torus of a large mass, rapidly accreting onto it. 
The most efficient cooling mechanism for the dense torus 
is provided by neutrino emission, and neutrino annihilation 
seems to be an important process for the fireball generation 
(even though it occurs on timescales that are typically too 
short to provide the bulk of the GRB energy; e.g. Ruffert 
& Janka 1999; 2001 but see also Setiawan, Ruffert & Janka 
2004). 
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The durations, T, of GRBs range from milliseconds to 
over a thousand of seconds, and are distributed in two dis- 
tinct peaks that constitute two GRB classes: short (T < 2 
sec) and long bursts (T > 2 sec; Kouvelietou et al. 1993). 
For the long bursts, signatures of an accompanying super- 
nova explosion have been detected in the afterglow spectra 
(Stanek et al. 2003). This strongly favors the "coUapsar" 
interpretation for their origin. Furthermore, the GRB posi- 
tions inferred from the afterglow observations are consistent 
with the GRBs being associated with the star forming re- 
gions in their host galaxies. However, all the observations 
that we have of the GRBs counterparts, i.e. afterglows and 
host galaxies, have been obtained for the long bursts (see the 
review by Zhang & Meszaros 2003) and no such signatures 
have been observed in the case of short bursts. Therefore 
the merger model is still thought to be viable for the case of 
short bursts, the nature of which remains mysterious. 

In the merger scenario, the duration of the GRB event 
is much longer than the dynamical timescale but rather 
comparable to the viscous timescale of the accretion disc. 
Therefore the phase of accretion sets the burst duration, 
differently from the coUapsar model where there is an ex- 
ternal reservoir of stellar matter which feeds the accre- 
tion torus. In general, accretion discs in the context of 
GRBs are expected to have typical densities of the order of 
10^""^^ g cm"^ and temperatures of lO" K within 10 - 20 
Schwarzchild radii (Rs = 2GM/(?). Thus, the accretion pro- 
ceeds with rates of a solar to several solar masses per second. 
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In this "hyper-accreting" regime, photons become trapped 
and are not efficient at coohng the disc. Neutrinos, however, 
are produced by weak interactions in the very dense and 
hot plasma ("neutrino-dominated accretion flow", NDAF, 
Popham Woosley & Fryer 1999; Narayan, Piran & Kumar 
2001; Kohri & Mineshige 2002; Di Matteo, Perna & Narayan 
2002). 

In this paper we compute the time-dependent evolu- 
tion of a massive accretion disc. We perform our calculations 
both within the context of a compact object merger scenario 
and in the coUapsar model, where the external supply of a 
matter, e.g. from the fallback from a supernova explosion 
(MacFayden, Woosley & Heger 2001), is essential to extend 
the burst duration into the range appropriate for long bursts. 
We investigate the evolution of the disc dominated by neu- 
trino cooling, and study the changes of neutrino and photon 
luminosities with time. In particular, we focus our discus- 
sion with the long term evolution (few tens of hours) of these 
discs and its implication for emission line models reported in 
the early X-ray afterglow of several GRBs. Our work is thus 
complementary to the 3D (or 2D) time-dependent numerical 
studies of Setiawan et al. (2004) and Lee & Ramirez-Ruiz 
(2002) in being able to extend to the study of the evolution 
of the bursts (in both coUapsar and merging scenario) to 
much longer timescale. The paper is organized as follows. In 
Section |^ we describe the basic assumptions of the model 
of hyper-Eddington accretion disc and the method used in 
time-dependent numerical simulations. In Sectionj^we show 
the evolution of the disc density and temperature, as well 
as the resulting neutrino luminosity. In Section 4 we discuss 
the implications of the long term evolution of the accretion 
disc for models of line production in GRB afterglow spectra. 
Finally, we summarize our results in Section |^ 



2 METHOD 

The structure of a steady-state, neutrino-dominated accre- 
tion flows (NDAFs) has recently been studied by Popham, 
Woosley & Fryer (1999), Kohri & Mineshige (2002) and Di 
Matteo, Perna & Narayan (2002). Here we use a similar pre- 
scription to calculate the initial disc configuration, which 
afterwards is subject to viscous evolution. 



2.1 Assumptions and model parameters 

Throughout our calculations we use the vertically integrated 
equations of the disc structure, however we note, that in 
most of the cases the disc becomes moderately geometri- 
cally thick {H ~ 0.5r). We assume that the disc extends not 
further out from the central black hole than 50 Rs, as shown 
e.g. in the recent simulations performed by Rosswog, Speith 
& Wynn (2004) . If the disc size were larger, the fiow would 
become highly convection-dominated instead of forming an 
NDAF (Narayan, Piran & Kumar 2001). The initial mass 
of such a disc is Mdisc ~ 47r_RoutEout, which is ^ 0.35Mq 
for the accretion rate M = 1 Mq/s. The angular veloc- 
ity of the disc is assumed to be Keplerian, Q = ^ GM j r'^ ., 

and the sound speed is Cs = \/ Pj p = n_ff. A non-rotating, 
Schwarzschild black hole is assumed and the inner radius of 
the disc is always at 3 i?s. 

For the disc heating we use the standard a viscosity 



prescription of Shakura & Sunyaev (1973) and we adopt a 
canonical value of a = 0.1. The total pressure P consists 
of the gas and radiation pressure as well as the pressure of 
degenerate electrons: 
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where we take the mass fraction of free nucleons Xnuc = 
30.97(p/10")(-^'"'Hr/10'°)(^''**) exp(-6.096/r/10i") if 

-'fnuc < 1.0 and Xnuc = 1.0 elsewhere, and we assume the 
molecular weight per electron /Xe = 2. The radiation pres- 
sure includes the contribution from the electron-positron 
pairs via the coefficient 11/12. The radiation pressure is 
always taken into account in the hydrostatic balance and 
thermodynamic relations, however we neglect it in the 
viscous heating. Therefore we have: 
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The cooling in the disc is due to advection, radiation 
and neutrino emission. The advective cooling in a stationary 
disc is determined from the global ratio of the total advected 
flux to the total viscously generated flux (e.g. Paczyiiski 
& Bisnovatyi-Kogan 1981; Muchotrzeb & Paczynski 1982; 
Abramowicz et al. 1988): 
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Here /3 is the ratio of the gas plus degeneracy pressure to the 
total pressure (3 — (-Pgas + Pdcg)/P- In the initial stationary 
disc we assume that gadv is approximately constant and of 
order of unity, but in the subsequent time-dependent evo- 
lution the advection term is calculated more carefully, with 
appropriate radial derivatives. 

In case of photon and electron-positron pairs in the 
plasma the radiative cooling is equal to: 

_3^ad£_iiar^ ,.. 

Orad - ^^ - -j;^ (8) 

where we adopt the Rosseland-mean opacity k = 0.4-f 0.64x 

lO^V^"^ cmVg. 

The most significant neutrino emission processes are: 
gNc = 9 X 103^(p/10"')(r/10")^Xnuc ergs/cm^s, the cool- 
ing rate in the non-degenerate URCA process (Qian & 
Woosley 1996); g^+o- = 4.8 x 10^^(r/10")^ ergs/cm^s, 
the emission due to the electron-positron pair annihilation 
(Itoh et al. 1989); gtroms = 1.5 x W'^'^ {T/10^^f"^{p/10^°f 
ergs/cm'' /s is the the bremsstrahlung emission in the non- 
degeneracy regime of nucleons (Hannestad & Raffelt 1998) 
and gpiasmon = 1.5 X 10^2(r/10")V(7p) ergs/cmVs is the 
plasmon decay rate (Ruffert, Janka & Schafer 1996). The 
largest contribution to the neutrino emission is due to the 
electron capture on nucleons gNo. The electron-positron an- 
nihilation rate could be neglected in the case of complete 
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degeneracy of electrons, i.e. for X^/kT >> 1, where Ac is the 
chemical potential of electrons. However, in the subsequent 
calculations the degeneracy parameter is never extremely 
large, and therefore gj.+o- is also taken into account. 

In the above processes either electron flavor neutri- 
nos and antineutrinos are produced (URCA process, plas- 
mon decay), or all the neutrino flavors, including the 
heavy lepton neutrinos and antineutrinos (pair annihila- 
tion, bremsstrahlung) . In this work we further neglect the 
possible anisotropics between neutrinos and antineutrinos 
as well as between different flavors, that could arise e.g. 
from a multiflavour neutrino leakage scheme (e.g. Rosswog 
& Liebendorfer 2003). 

We treat, in a simplified scheme, neutrino cooling both 
in the optically thin and optically thick regimes. 

In the (i) optically thin regime the cooling is simply given 
by: 
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In the (ii) optically thick regime we include in our calcu- 
lations the prescription for the neutrino opacity according to 
Di Matteo et al. (2002). This will allow us to switch smoothly 
between the optically thick and thin regimes. Each neutrino 
emission process described above has an inverse process cor- 
responding to absorption. Absorption onto protons or onto 
neutrons (the inverse of the URCA process) gives rise to the 
large majority of the (a) absorptive optical depth, thus given 
by (see Di Matteo et al. 2002 for details): 

QNeH 
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A very important contribution to the overall opacity for 
all neutrino species comes from neutral-scattering off nucle- 
ons. The (b) scattering optical depth is given by (see Di Mat- 
teo et al. 2002 for derivation of scattering optical depth): 

Tsc = SpusH = 8.1 X W~''{T/10^^f{p/W^°)H (11) 

where the factor 3 assumes equal contribution for all neu- 
trino species, i.e. thermal equilibrium. 

The neutrino cooling in this case is then given by: 

O; = — ^^ r- (12) 
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In addition, in this regime, the equation of state should con- 
tain the neutrino pressure: 
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Finally, the entropy density due to neutrinos. Si, = 4/3 x 
7/8ar'', is included in the advective cooling when we calcu- 
late the initial disc configuration. 

The last equation that is essential to close our set of 
equations for the stationary disc is the total energy flux 
-Ftot- It is dissipated within the disc at a radius r and is 
determined by the global parameters: 



-Ftot — 
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(Note, that the above equation will not be used in the subse- 
quent time-dependent calculations.) In order to calculate the 
initial stationary configuration, we solve the energy balance: 



-Ftot = Qtisc = Qldv + Orad + Q^- In this energy balance 
we did not include the cooling term due to the photodis- 
integration (Di Matteo et al. 2002), while we include the 
radiative cooling. We checked that this term is much less 
than the neutrino cooling rate in the innermost parts of the 
disc, when the hot disc is cooled by neutrinos, and less than 
advective and radiative cooling, when these two mechanisms 
start playing important role. 

The initial configuration of the disc is calculated by 
means of a simple Newtonian method. This lets us deter- 
mine the radial profiles of density and temperature, as well 
as the disc thickness and opacity at time t = 0. 



2.2 Time evolution 

In the previous subsection we described the initial, station- 
ary disc configuration. Having computed the initial disc 
state we allow the density and temperature to vary with 
time. We solve the time-dependent equations of mass and 
angular momentum conservation for such a disc: 
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and the energy equation: 
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Here E = Hp is the surface density and v,- ~ (3i^)/(2r) is 
the radial velocity in the disc, i/ = {2Pa)/{3pQ,) is the kine- 
matic viscosity. The cooling term Q~ consists of radiative 
and neutrino cooling, given by Equations 151 and lUl for IT^ . 
The advection is included in the energy equation via the ra- 
dial derivatives. In the heating term Q+ we include gas and 
degeneracy pressure (in the optically thin regime) and also 
neutrino pressure (in the optically thick regime). Note, that 
this equation includes the radiation entropy, via the coeffi- 
cient /3, which is not equal to unity due to the presence of 
radiation pressure in the total pressure. 

We solve the above set of time-dependent equations us- 
ing the convenient change of variables, y — 2r^'^ and H = 
j/E, at the fixed radial grid, equally spaced in y (see Janiuk et 
al. 2002 and references therein) . The number of radial zones 
is set to 20. After determining the solutions for the first 
100 time steps by the fourth-order Runge-Kutta method, 
we use the Adams-Moulton predictor-corrector method, al- 
lowing the time-step to vary, when needed. 

We choose the no-torque inner boundary condition, 
Ein = Tin = (see Abramowicz & Kato 1989). The outer 
boundary of the disc is parameterized by an external accre- 
tion rate A/ext in case of the collapsar scenario. 

In case of the merging compact objects scenario there 
is no matter supply to the disc, and we can have Eout = 
Tout ~ 0. However, the outer edge of the disc should ex- 
pand to conserve the angular momentum, and in principle 
the density would vanish at the infinite distance. In order 
to satisfy this requirement, we supplement our radial grid 
with 10 empty rings between 50 and 96 -Rs, while the outer 
boundary condition E = T = is located at 96 -Rs- These 
empty rings gradually fill with material that sinks from the 
outermost disc radius by means of the viscous diffusion, and 
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take the angular momentum form the disc. The extension of 
the empty zone is large enough to ensure that at the end of 
the calculations the surface density in N-1 ring is less than 
5% of the surface density in the zone 20 (the edge of the 
initial disc). 



3 RESULTS 

We first analyze the local evolution of the disc on the T — E 
plane, and then we show examples of time evolution of the 
neutrino luminosity for a chosen set of parameters. 



3.1 Stability curves 

A few examples of stability curves calculated at 4, 10, 20 
and 46 Rs for the mass of the central object equal 3Mq are 
shown in Figure^ The open points represent the stationary 
solutions, calculated from the grid of steady-state models 
(optically thin regime). The solid points represent the sub- 
sequent solutions of the time-dependent model. Deviation 
between the stationary and time-dependent solutions arises 
from the simplified description of advection, used in the sta- 
tionary model. In that model we had a constant value of the 
advection parameter given by Eq. |7| qadv = 1.0 and does 
not depend on radius, while in the time dependent model 
the radial derivatives are calculated properly. However, this 
simplification of the stationary model does not influence the 
time-dependent results, since we start the evolution from 
the upper, neutrino-cooled branch, where advection term is 
not crucial. 

We see that the resulting curves display the pattern 
already discussed in Kohri & Mineshige (2002). Both up- 
per and lower branches of different slopes are viscously and 
thermally stable. On the upper branch the dominant cooling 
mechanism is the neutrino cooling, and the radiation pres- 
sure is much lower than gas and degeneracy pressures (most 
of this branch is in fact gas-pressure dominant, and the de- 
generacy pressure can exceed it only for very large densities, 
p > 3 X lO^'^ g/cm^). On the lower branch the radiation 
pressure dominates, however the flow is cooled by advection 
rather than radiation, and therefore the disc should remain 
stable. Below this branch, for much lower densities and tem- 
peratures, we expect the unstable (radiatively cooled and 
radiation pressure dominant) branch to appear, as shown in 
detail by Pringle, Rees & Pacholczyk (1974) and Lightman 
&Eardley (1974). 

We start the disc evolution from the top of the upper 
branch, and first we let it cool down very quickly due to 
the zero mass inflow at the outer edge. In the initial con- 
figuration, neutrino cooling is dominant in the entire disc. 
Although the outer boundary condition imposes a fast den- 
sity decrease at the outer edge of the disc, the inner parts 
do not respond immediately to that. For the first ~ 1 sec- 
ond the disc annuli evolve slowly along the upper branch, 
and the density of the disc decreases as the matter falls into 
the black hole. In this branch the temperature changes are 
not very significant, so it is the substantial density evolution 
that drives the steep fall of the disc luminosity. 

The evolution starts to proceed faster when all the disc 
annuli find themselves on the lower branch, where advective 
cooling becomes dominant. As the whole disc is dominated 
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Figure 1. Local evolution of the disc on the surface density - 
temperature plane, plotted for 4 values of radius: 4.06, 9.85, 20.64 
and 46.1 Rs, in case of the disc optically thin to neutrinos. The 
open triangles mark the stability curve resulting from stationary 
disc solutions. The solid points are the time-dependent solutions. 



by advective cooling, we reach the so-called slim disc branch. 
Here the radiation pressure is larger than gas and degener- 
acy pressure. As a result of this, we may expect the outer 
disc edge to develop an instability leading to the disc break 
down, if the viscous heating was assumed proportional to 
the total (gas and radiation) pressure. On the other hand, if 
the viscous heating is assumed proportional only to the gas 
pressure, there is no problem with instabilities during the 
slim disc evoluton. Therefore we choose this prescription, 
also having in mind the recent results of 3-D simulations 
performedd by Turner (2004) , which show that the radiation 
pressure instability is strongly suppressed when the angular 
momentum transport in the accretion disc is described by 
MHD calculations, without the a parameter. This heating 
prescription makes no difference at the NDAF branch, where 
more important contribution is given by gas, degeneracy and 
neutrino pressure, which we always take into account. Since 
neutrino cooling is no longer important (L^ < 10^" ergs/s), 
further calculations would not be relevant for the the study 
of GRB central engine and we stop our calculations here. 

The disc evolution proceeds much slower when the ac- 
cretion rate at the outer edge is non-zero. The phase of the 
most efficient neutrino cooling, i.e. when the innermost disc 
parts are cooled mainly by the neutrino emission and not 
by advection, is extended to ~ 10 seconds. Also the sec- 
ond phase, of the advection-dominated disc, is much longer. 
The disc evolves along the advective branch more slowly and 
there is no rapid emptying of the outer parts. 
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Figure 2. The time evolution of the neutrino luminosity in case 
of no external mass supply to the disc (optically thin case). The 
initial accretion rate is A/ = 1 A/q/s. 



3.2 Lightcurves 

The outer boundary condition for the disc evolution depends 
on the accretion rate, Mcxt- If we assume that the torus was 
formed during the merger of two compact objects, NS-NS 
or BH-NS, it will consist of the debris matter from the dis- 
rupted neutron star and there will be no source of additional 
material infalling onto the disc. On the other hand, if we take 
into account the "coUapsar" model, the associated supernova 
explosion should provide a declining fallback of matter, that 
accumulates mostly in the equatorial plane. The variations 
of the mass supply to the accretion flow will strongly in- 
fluence the disc evolution and therefore result in significant 
changes in the lightcurves. 

The lightcurves represent the neutrino luminosity of the 
disc. The luminosity is given by: 



L. 



Qy 2TTrdr 



(17) 



where Q is given by either Eguationl^or H^ as appropriate 
for the optically-thin and thick regimes. 



3.2.1 Optically thin regime 

In Figure |21 we show an example of a neutrino lightcurve 
calculated under the assumption of an optically-thin neu- 
trino disc for the case of zero mass supply. Other parameters 
were always the same: mass of the central object M = 3Mq , 
a = 0.1, -Rout = 50-Rs, and the starting accretion rate was 
Mo = 1 Mq/s. 

In Figure 121 we show the neutrino lightcurve calculated 
for the case of additional matter supply at the outer edge of 
the disc, which declines with time as: 



Figure 3. The time evolution of the neutrino luminosity in case 
of additional external mass supply to the disc, declining with time 
as in the canonical supernova fallback prescription, Afoxt ~ t~^' ^ . 
The initial accretion rate is Af = 1 Mq/s. The disc is optically 
thin to neutrinos. 



Mext = Mot-"^^ 



(18) 



This is the time dependence of the fallback of material de- 
rived for the supernova SN 1987A (Chevalier 1989) and we 
adopt it as a canonical rate of the fallback in case of GRBs 
associated with supernovae. The starting accretion rate was 
again Mo = 1 Mq/s. 

Clearly, the supply of material at the outer edge of 
the disc extends the duration of the hot phase, when the 
neutrino cooling dominates. GRBs associated with a super- 
nova explosion can therefore have much longer durations, 
of the order of tens of seconds, depending on the initial 
mass of the disc and accretion rate. The neutrino lumi- 
nosity declines exponentially, with some slight luminosity 
fluctuations in the initial phase, when the neutrino cool- 
ing remains dominant at least in the innermost disc parts. 
The lightcurve smoothens as soon as the whole disc switches 
to the advection-dominated mode, but its slope does not 
change much. 

3.2.2 Optically thick regime 

When the disc is optically thin to its neutrino emission, the 
cooling rate due to neutrinos, Q^, is given by Equation [51 
However, the disc can become optically thick to neutrinos, 
especially in its innermost regions and for high accretion 
rates. The optical depth is larger than unity in the initial 
disc state. Although it drops very fast in the outer disc radii, 
still remains quite high for most of the disc evolution in 
the innermost parts. In this case the cooling rate due to 
neutrinos Q~ is given by Equation 1121 

In Figure |1] we show the neutrino luminosity resulting 
from the model in which we include the treatment of neu- 
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Figure 4. The time evolution of the neutrino luminosity in case 
of no external mass supply to the disc and neutrino optical depth 
included in the model. The initial accretion rate is M = 1 Mq/s. 



Figure 5. The time evolution of the neutrino luminosity in case 
of no external mass supply to the disc (compact objects merging 
scenario). The solid line shows the model with neutrino optical 
depth, while the dashed line is the optically thin case. The initial 
accretion rate is M = 1 A/q/s and a = 0.1. 



trino optical depths and add the neutrino pressure term in 
the equations. As expected, the most significant change in 
comparison with the optically thin case (Figure 2) occurs at 
the very beginning of the evolution, when the innermost disc 
is optically thick to its neutrino emission. After about 0.2 
seconds most of the disc becomes optically thin to neutrinos 
and the solution is almost the same as in Section f3.2.1l 

Recently, Lee et al. (2004) computed the dynamical 
simulation of a short GRB, resulting from the merger of 
two neutron stars. They showed that inclusion of the neu- 
trino optical depth influences the duration and variability 
of the short burst. Our simulation confirms this conclusion: 
in Figure |K| we show the comparison of the neutrino lumi- 
nosity evolution in the optically thin and thick regimes, in 
the zoom-in emphasizing the shortest timescale by means of 
the logarithmic units. As the figure shows, in the optically 
thick case the neutrino luminosity around 20 milliseconds 
is enhanced, consistently with Lee et al. simulations. Then, 
the luminosity declines slightly steeper than for the thin disc 
and the duration of the burst is shorter. 

In Figure |S] we plot the disc evolution on the surface 
density vs. temperature plane in case of the disc optically 
thick to neutrinos. The main difference in comparison to 
the optically thin case (Figure^ is at the inner radii, where 
the local solutions achieve higher temperatures, T > 10^^ K. 
Further out, and also later in time, when the disc becomes 
advective and neutrino optical depth is less than 1.0, the 
solutions match the optically thin case. 

In case of a long burst, the difference between the opti- 
cally thin and optically thick regimes is not significant, since 
after ~ 1 second of the disc evolution the both optical depths 
Tabs and Tsc become less than 1.0 and we can substitute the 
optically thick equations with the optically thin ones. Thus, 
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Figure 6. Local evolution of the disc on the surface density - tem- 
perature plane, plotted for 4 values of radius: 4.06, 9.85, 20.64 and 
46.1 Rg, in case of the neutrino-thick inner disc. The points are 
the time-dependent solutions, for the time range - 0.4 seconds. 
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in practice, the resulting lightcurve does not differ from tfiat 
shown in Figure |21 



4 LONG-TERM EVOLUTION OF THE 
ACCRETION DISC AND ITS 
IMPLICATIONS FOR EMISSION LINE 
MODELS 

Our time-dependent code allows us to follow the evolution 
of the GRB accretion disc up to times much longer than the 
duration of the GRB itself. This is valuable information in 
light of the peculiarities observed in the early stages of some 
afterglows (e.g. Fox et al. 2003), which hint at a more ex- 
tended period of energy injection than the GRB phase itself. 
What is however particularly relevant for, is the interpreta- 
tion of the narrow emission features that have been reported 
in the early X-ray afterglow of several bursts (Yoshida et al. 
1999; Antonelh et al. 2000; Piro et al. 1999, 2000; Reeves 
et al. 2002). Although the statistical significance of the de- 
tections is not compelling, the interpretation of these lines 
bears important implications for our understanding of GRB 
progenitors. 

The emission lines are detected at the rest frame fre- 
quency of the host galaxy, and therefore they cannot be 
produced within the fireball, which, during the X-ray after- 
glow phase, is still moving towards the observer at a very 
high speed, with a Lorentz factor F > 10. Any emission 
model for the lines must be confronted (as a minimum re- 
quirement) with the energetics of the lines, which have a 
luminosity (for isotropic emission) ~ lO*** — 10*^ erg/s and 
last for several hours. This luminosity values must be aug- 
mented by a factor I/77, where rj is the efficiency for the 
reprocessing of the impinging flux into line photons. The ef- 
ficiency of conversion of the X-ray ionizing continuum into 
Ka line photons was investigated by Lazzati et al. (2002). 
They assumed that the line is produced by refiection off 
an optically thick slab of material, illuminated by a power- 
law continuum. Under these conditions, which provide the 
most efficient way of reprocessing continuum photons into 
Ka lines, they found that the efficiency rj can be at most 
2%. This implies that the impinging continuum flux must 
carry an energy > 10^^ — 10^^ erg/s. 

This energy requirement could in principle be allevi- 
ated if the line emission were coUimated. In the model of 
Rees & Meszaros (2000), where lines are produced by mul- 
tiple scattering off the walls of an evacuated funnel along 
the rotation axis of the collapsing star, there can indeed 
be some coUimation of the emission, enhancing its observed 
flux. The importance of this effect was however shown to be 
not too significant by Ghisellini et al. (2002). In particular, 
they found that, if the main scattering opacity within the 
funnel is provided by free electrons, then photons scattered 
deeper in the funnel will be lost, and the resulting line am- 
plification can be at most a factor of 2. On the other hand, if 
the line under consideration is a resonant line, these photons 
may be rescattered by ions in the funnel wall and redirected 
towards the open space. In this case, line amplification can 
get up to a factor of 10 for an opening angle of the funnel 
of ~ 5 deg (and it gets smaller for larger angles). Therefore, 
even under the most favorable circumstances, the required 
luminosity of the continuum has to be > 10^^ — lO''^ erg/s. 
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Figure 7. Long-term evolution of photon (solid line) and neu- 
trino (dashed line) luminosities of the remnant torus, fed by the 
material fallback at its outer edge as in Eq. 1181 



The issue therefore remains what is that provides the energy 
to power the lines. 

The recent detection of a supernova associated with 
GRB 030329 (Stanek et al. 2003; Hjorth et al. 2003) gave 
strong support to the ideas that GRBs are the result of ex- 
plosion of massive stars. In the coUapsar model (MacFadyen 
& Woosley 1999), the remnant compact object is a black 
hole, and the "engine" is an accretion disc like the one we 
have studied here. The energy output at later times (i.e. 
after the GRB phase) is in this case given by the sluggish 
drain of orbiting matter into the newly formed black hole, 
after the supply of matter from the envelope of the collaps- 
ing star is exhausted. Our calculations allow us to evaluate 
the energy output from this remnant disc at the times (of 
the order of tens of hours) at which the emission features 
have been observed. Figure shows the evolution of the 
photon and neutrino luminosities from the very early times 
of the GRB phase to several days after the GRB. While 
the neutrino luminosity drops following the decrease in the 
accretion rate (Figure |SJl, the photon luminosity increases 
with time. This is due to the decreasing photon opacity of 
the disc, which overcomes the decrease in energy fiux that 
would be expected as a result of the declining mass accretion 
rate. The Rosseland mean opacity, calculated for the elec- 
tron scattering and free-free transitions (see Equation IHJ is 
very large in the initial disc, which is therefore almost totally 
opaque to photons (k > 3 x 10^ for r < 10-Rs). However, the 
opacity drops quickly and after ~ 10 minutes it approaches 
the value of 0.4. Therefore as long as the main source of 
cooling is the neutrino emission, the photon luminosity rises 
with time. The trend will reverse when L^ < L-y and the 
photon emission becomes the more important cooling mech- 
anism (note however that advection is still very large). The 
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Figure 8. Long-term evolution of tiie accretion rate in the torus, 
plotted for several radii: 5.27 -Rg (solid line), 9.85 Rg (dotted 
line), 20.64 Rs (short-dashed line), 35.37 Rs (long-dashed line) 
and 46.1 -Rs (dot-dashed line). 



decreasing mass accretion rate will now result in lower and 
lower photon emission from the disc. 

As the figure shows, at the times corresponding to the 
line observations, the photon luminosity is ~ 10''^ — 10^^ 
erg/ sec, several orders of magnitude below what is required 
to produce the emission features at the observed luminosity 
level. The simulations in Figures Q and |S] are performed for 
an initial accretion rate of 1 Mq/s, but they give similar 
results for any initial value of M that produce a GRB at 
early times. Since for higher initial accretion rates the neu- 
trino cooling of the disc is even more effective, the photon 
luminosity at late times is lower: for 2Mq/s it is 3.2 x 10^^ 
erg/s at f ~ 1 day, while for 1Mq/s it was 4.7 x 10''^ erg/s. 
The viscosity in the disc does not have a big impact on the 
long-term results either: the dissipation rate and the disc 
luminosity increase only slightly with a and the change is 
of the order of AlogL.^ = 0.3 and AlogL,^ — 0.5 when 
we change a form 0.1 to 0.3. Therefore, our simulations im- 
ply that the energy output from the accretion disc alone is 
not sufficient to power the GRB emission lines (if these are 
indeed real). 

An alternative possibility (Rees & Meszaros 2000) for 
producing a long-lasting energy output is by the collapse 
of a massive star into a fast millisecond pulsar (Usov 1994; 
Thompson 1994; Wheeler et al. 2000). In this case the late 
time energy is provided by the electromagnetic losses of the 
pulsar which, at a day after the burst, can easily be as high as 
10*'^ erg/s. Our results favor this scenario. Identification of 
the remnant compact object after the GRB explosion would 
provide an independent test of this scenario, and this should 
be possible in local galaxies after identification of GRB rem- 
nants (Loeb & Perna 1998; Efremov et al. 1998; Perna et al. 
2000). 



5 DISCUSSION 

The bimodal distribution of durations of GRBs is most prob- 
ably due to the different origin of the short and long bursts, 
and the proposed models of a GRB progenitor include bi- 
nary mergers or collapse of a massive star. In both cases 
the formation of a dense accretion disc around a black hole 
is expected, and one of the possible ways of extracting the 
energy from the central engine is the neutrino-antineutrino 
annihilation. This process can give rise to the relativistic 
fireball, if the neutrinos convert into electron-positron pairs 
in the region of relatively low baryon density. 

The merger scenario seems to be appropriate only in 
the case of short GRBs, of duration < 2 seconds. Due to the 
need for a low baryon density above the poles, the favored 
case is black hole - neutron star merger. The neutron star 
can either be catastrophically disrupted in the first approach 
to the black hole or transfer mass to the BH 2 or 3 times, 
depending on the initial mass ratio. In either case in the end 
it forms a torus, whose mass is between 0.25 and 0.7 Mq, 
and accretion rate of several AIq/s. Therefore the accretion 
time is of order of 0.1 - 1.5 sec (Janka et al. 1999). 

In the coUapsar scenario the black hole formation may 
be accompanied by a fallback of material, which over a pe- 
riod of minutes to hours feeds the accretion disc. If the initial 
accretion rate is sufficiently high to establish a dense disc 
cooled by neutrinos (M Js 0.1 Mq/s), the powerful, long 
GRB may be produced. 

The neutrino annihilation is most efficient near the ro- 
tation axis (Jaroszyriski 1993). This process is able to ac- 
count for the energy released in GRBs if a moderate beam- 
ing, f7/27r ~ 10~^, is involved (Ruffert & Janka 1999). In 
this case the isotropised energy of a burst is of the order 
of 10^^-10^^ erg/s. The jet coUimation mechanism is pos- 
sible e.g. due to the neutrino-driven wind, which can be 
expected from the outer disc regions if the neutrino lumi- 
nosity is Li, > lO^'^ ergs/s and the estimated wind power 
is Lw = 2 X 10*^(L,./10^^)(^'^''^) (Rosswog & Ramirez-Ruiz 
2003). 

In addition to the fraction of accretion energy that goes 
into the jet via neutrinos, the energy extraction from the 
central engine is possible via magnetic fields (Blandford & 
Znajek 1977; Meszaros & Rees 1997; Kluzniak & Ruderman 
1998). In this case, isotropic luminosities as high as 10^^- 
10^"^ erg/s are possible, even without beaming. The burst 
energy can also be increased due to the effects associated 
with the Kerr metric of a rotating black hole. 

Here we start our calculations when the accretion disc 
hass already established its quasi-steady state and evolves 
in the viscous timescale. We neglect the first, very violent 
stage of the disc formation, during which the material within 
a free fall time settles on the equatorial plane. The centrifu- 
gal force is balanced by gravity and the torus forms with a 
surrounding accretion shock. Inside this shock the temper- 
atures and densities are high enough to provide anhanced 
neutrino emission during a few seconds before the disc for- 
mation (MacFadyen & Woosley 1999). 

We computed the time dependent evolution of the ac- 
cretion disc believed to be the GRB central engine. Our 
simulations took into account a physical equation of state, 
including gas and radiation pressure as well as the pressure 
from degenerate electrons. We calculated the energy balance 
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between the viscous heating and neutrino losses due to ther- 
mal processes and neutronization, in both optically thin and 
thick regimes. We also took into account the advective and 
radiative cooling. Time-dependent 2D hydrodynamical sim- 
ulations of an accretion disc in GRBs have recently been 
calculated by Lee & Ramirez- Ruiz (2002), Lee et al. (2004) 
and Setiawan et al. (2004). Here we use the vertically inte- 
grated equations and the simulations are one-dimensional, 
which on one hand results in a great simplification of the 
problem, but on the other hand allows us to get insight into 
the various physical processes that act during the disc evo- 
lution. The most valuable aspect of such an approach is the 
possibility to study the long-term evolution of the remnant 
disc, especially in the frame of a long burst (coUapsar) sce- 



The short bursts are powered by the neutrino cooled 
accretion disc, which exists for about ~ 2 seconds. At this 
phase the neutrino luminosity remains very high during 
the first ~ 1 — 1.5 seconds, while dropping as soon as the 
whole disc becomes advection dominated. In case of the long 
bursts, the fallback of material that occurs during the ac- 
companying supernova explosion plays an essential role. The 
neutrino luminosity remains very high during the first 2-3 
seconds of the burst, and then drops exponentially, being 
sufficiently high to feed the GRB during several tens of sec- 
onds. 

The hyper-accreting, neutrino cooled disc is in either 
case very opaque to photons and the photon luminosity is 
always negligible whenever the neutrino emission is signif- 
icant. However, as the disc evolves and the neutrino cool- 
ing becomes less efficient, the disc opacity due to electron 
scattering and free-free transitions decreases, thus allowing 
more and more photons to escape. Therefore the photon lu- 
minosity rises in time up to the level at which it exceeds 
the decreasing neutrino luminosity. Later on, since the disc 
is cooled mainly by photons (and advection), and during its 
evolution the accretion rate and energy generation decrease, 
the photon luminosity drops with time. 

This stage of significant photon emission from the disc is 
extremely short and therefore not very important in case of 
the binary merger model. However, in the coUapsar scenario 
it is much more extended in time and can last over several 
days to weeks, depending on the initial accretion rate and 
viscosity (the greater initial M and a, the longer the photon 
emission phase) . Our calculations allow us to determine the 
residual level of photon luminosity from the disc at the times 
(~ tens of hours) corresponding to the claimed observations 
of X-ray emission features in a few GRBs. 

The emission from the disc could in principle power 
the emission lines, if the luminosity were sufficiently high. 
However, in our model this would require quite low initial 
accretion rates and an extremely high value of the viscos- 
ity parameter in the disc (for a > 0.9 and A/ = O.IMq/s 
we could expect the X-ray luminosity of the order of 10** 
erg/s), whereas for a standard value of viscosity a — 0.1, 
which is commonly adopted, the disc alone is not capable 
to power the lines. Therefore our simulations would rather 
favor the possibility of a long-lasting energy release by the 
electromagnetic losses of a newly born millisecond pulsar. 
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